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High- Resolution Measurements Of Intersystem Bands Of Carbon 

Monoxide Toward X Persei 

Yaron Sheffer 1 , S. R. Federman 1 , and David L. Lambert 2 

ABSTRACT 

In an echelle spectrum of X Per acquired with the Space Telescope Imaging 
Spectrograph we have identified individual rotational lines of 11 triplet-singlet 
(intersystem) absorption bands of 12 CO. Four bands provide first detections for 
interstellar clouds. From a comparison with the £ Oph sight line we find that X 
Per is obscured by a higher 12 CO column density of 1.4 x 10 16 cm -2 . Together 
with the high spectral resolution of 1.3 km s -1 , this allows (i) an improved mea- 
surement of previously published /-values for seven bands, and (ii) an extraction 
of the first astrophysical oscillator strengths for d— X (8—0), (9—0), and (10—0), 
as well as for e—X (12—0). The 13 CO d—X (12—0) band, previously suspected 
to exist toward £ Oph, is now readily resolved and modeled. Our derived inter- 
system /-values for 12 CO include a few mild (< 34%) disagreements with recent 
predictions from a perturbation analysis calculated for the interstellar excita- 
tion temperature. Overall, the comparison confirms the superiority of employing 
multiple singlet levels in the calculations of mixing coefficients over previous 
single-level predictions. 

Subject headings: ISM: abundances - - ISM: molecules - - molecular data - 
ultraviolet: ISM 



1. INTRODUCTION 

Observers of interstellar sight lines are familiar with the strong A 1 II— X X S + permit- 
ted absorption bands of CO below 1545 A. Less familiar to observers are the intersystem 
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bands of CO, which involve the triplet states a' 3 E + , d 3 A;, and e 3 £~, whose electric dipole 
transitions to the ground state X 1 Y, + are spin-forbidden. These triplet states have poten- 
tial energy curves that cross the curve of A 1 I1 — see Figure 1 of Rostas et al. (2000) — and 
via mutual perturbations acquire a small percentage of 1 I1 character at the expense of the 
singlet state. The resulting triplet-singlet /-values are relatively small, thus turning the 
intersystem bands into very useful astrophysical probes along interstellar (IS) sight lines for 
which the A—X bands are optically thick, see Morton & Noreau (1994, hereafter MN), Fed- 
erman et al. (1994, hereafter F94). To exploit such probes, accurate /-values of intersystem 
transitions are required. Rostas et al. (2000, hereafter Rostas et al.) presented predictions 
from extended calculations of the triplet-singlet mixing coefficients, partly as an attempt to 
resolve discrepancies between previous theoretical (MN) and astrophysical /-values. Rostas 
et al. also conducted low-resolution laboratory measurements at room temperature to im- 
prove upon the experimental results of Eidelsberg et al. (1992). High-resolution laboratory 
measurements on two intersystem bands were also reported by Stark et al. (2002). 

Using the Hubble Space Telescope^ (HST), we obtained a high- resolution echelle spec- 
trum of the 09.5V star X Per (HD 24534) with the Space Telescope Imaging Spectrograph 
(STIS). Permitted A—X bands along this line of sight were first modeled by Kaczmarczyk 
(2000a), who based his analysis of 12 CO on archival HST Goddard High- Resolution Spec- 
trograph (GHRS) observations. For our STIS observations, the weakest permitted band 
available for analysis is the A—X (8—0) band, which still has a large optical depth (r) in 
the cores of the strongest lines, t(R(0)) ~ 70. Therefore, the column density (N) of 12 CO 
toward X Per can be measured more reliably from the observed intersystem bands, which 
have a much smaller r. In this Letter, we determine iV( 12 CO) through a comparison with 
our previous observations of intersystem bands toward £ Oph (F94). Then we derive IS 
/-values toward X Per for all seven "old" bands and for five intersystem bands seen in IS 
spectra for the first time. Finally, we compare our /-values with the predictions of Rostas 
et al. Regarding nomenclature, we make use in this paper of the short notation introduced 
by Rostas et al., where, e.g., d!2 stands for d—X (12—0). 



3 Based on observations obtained with the NASA/ESA Hubble Space Telescope through the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS5-26555. 
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2. OBSERVATIONS AND MODELING 

STIS spectra of X Per were acquired in 2001 February and March with grating E140H 
during 10,728 s over five orbits, yielding the data sets o64812010-030 and o64813010-020. 
The star was observed through the smallest aperture (0.1" x 0.03"), which was designed 
to provide the highest resolving power, R = 200,000 (Jenkins & Tripp 2001). Elaborate 
reduction procedures had to be employed in order to correct for odd-even sub-pixel effects 
in the unbinned (2048-pixel) exposures. Normal reduction routines were then followed by 
processing the spectra in the STSDAS environment of IRAF. The five orbital exposures were 
co-added into a single spectrum of 42 orders, sporting a signal-to-noise ratio of 70, or 100 for 
features detected and combined from two adjacent overlapping orders. Since the continuum 
is modulated by shallow photospheric features from X Per, short spectral segments centered 
on features (bands) of interest were rectified and normalized before being compared to a fitted 
synthetic model of the pertinent feature. Flux and wavelength calibrations were performed 
by the Space Telescope Science Institute pipeline and supplied with the binned, 1024-pixel 
exposures. 

A previous analysis of 12 CO by Kaczmarczyk (2000b) inferred that this line of sight 
involves two unresolved velocity components. However, the significantly higher resolution 
of our spectrum shows profile asymmetries due to underlying structure, which required four 
components for a good fit. This is the same number of components found in preliminary 
fits of high-resolution K I spectra of X Per (D. E. Welty, private communication), although 
we derived a different set of relative velocities, fractions, and component widths. Some 87% 
of the total iV( 12 CO) resides in a main component with a Doppler parameter (fr-value) of 
0.39 km s^ 1 . The other three components have between 2% and 9% shares of N, with 
6-values between 0.17 and 2.09 km s -1 , as will be fully described in a future paper. The 
highly saturated A— X bands could be fitted only once the total N was established via the 
intersystem bands. Our method of analysis employed a multi-parametric spectrum synthesis 
code based on Voigt profile line transfer equations given in Black & van Dishoeck (1988). 
The computed absorption spectrum was convolved with a Gaussian instrumental profile 
(Jenkins & Tripp 2001) and then matched with the data via rms-minimization down to 
relative parameter steps of 10~ 4 or less. 



3. THE COLUMN DENSITY OF 12 CO 

In order to determine /-values for intersystem bands, iV( 12 CO) toward X Per must be 
known. For optically thin features, the equivalent width {W\) is proportional to r, i.e., 
W\ oc fN. We can, therefore, utilize the published value of N toward ( Oph and the ratio 
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of W^s for optically thin bands toward the two stars in the determination of N toward 
X Per. The result should be independent of the intersystem /-values and of the modeled 
cloud structure, provided we exclude the W\ ratio for the optically thick a' 14. Column 4 of 
Table 1 lists the ratios of band W\ toward X Per over band W\ toward ( Oph. Among the 
optically thin bands, the highest W\ ratio belongs to e4, which was the faintest intersystem 
band analyzed by F94. In addition, e4 also sports the largest deviation of derived /-value 
in the low-r sample of F94. It is clear that F94 underestimated the low-resolution W\ for 
e4 due to placing the continuum too low (see their Figure 1). Dropping the e4 band, and 
thus avoiding a 5% increase in the average, we obtain W\(X Per)/W,\(C Oph) = 5.54 ± 
0.87. Under the necessary and fulfilled condition of optical thinness, this is also the ratio of 
iV(X Per) over N(( Oph). Here we shall adopt the Lambert et al. (1994) column density 
toward ( Oph, N( 12 CO) = (2.54 ± 0.16) x 10 15 cm' 2 , which was also employed by F94. 
Based on the measured N ratio, we infer that iV( 12 CO) toward X Per is (1.41 ± 0.22) x 10 16 
cm" 2 . Kaczmarczyk (2000b) found (1.0 ± 0.2) x 10 16 cm" 2 from an analysis of permitted 
A—X bands; the agreement between the two determinations is at the 2-a level. Recall, 
however, that Kaczmarczyk's model has fewer cloud components and larger 6-values, readily 
accounting for its reduced column density. 



4. BAND OSCILLATOR STRENGTHS 

The crossings of vibrational levels in a triplet state and vibrational levels in the A 1 U 
state occur near certain rotational J' levels. Therefore, mixing coefficients depend on J' and 
thus intersystem / band -values generated by interaction with A 1 I1 have a dependence on T ex . 
Whereas all previous calculations of mixing coefficients focused on the vibrational level of 
the A state (va) closest to the perturbing triplet level, i.e., a single-^ treatment, Rostas 
et al. were the first to use multi-iu computations that included all A levels with v' — 
to 12. They derived /band-values that for some intersystem bands were in better agreement 
with their own experimental results (T ex = 300 K), and/or with the IS results of F94 (T cx = 
4.2 K). Still, the calculated /-values of Rostas et al. are "a few percent" off, because the 
changes in f v i v ii of the "parent" A—X bands induced by mixing with other triplet levels were 
ignored. In addition, differences of a few percent are expected due to the adopted sources for 
A—X /-values: MN and we used values from Chan, Cooper, & Brion (1993), while Rostas 
et al.'s preference lies with values published by Eidelsberg et al. (1999). 

The MN line oscillator strength, fjr j», which we use as input for the fits, includes f v ' v " 
for the permitted A—X band and the usual Honl-London rotational factor, as well as the 
single-iu strength of the perturbation. According to Rostas et al., for triplet levels that 
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cross the A levels at low values of J', single-tu and multi-f a calculations should not differ 
by a large factor since the interaction with the nearest "parent" va level is predominant. 
Therefore, using fjijn -values from MN is a valid approach while searching for relatively 
small corrections due to multi-iu interactions. This condition is appropriate for all bands 
studied here, except for d7 , dlO, and to a lesser extent, e4. On the other hand, intersystem 
bands that are farther from their nearest parents have "J-independent" /band-values due to 
their high- J' crossing, especially under IS conditions. Again, using /j/j/ /-values from MN 
poses no difficulties for scaling the integrated band /-values thanks to negligible high- J" level 
populations. Indeed, we made no attempt to fit rotational lines independently of each other; 
our syntheses did not indicate a need to adjust individual /j/j/ /-values. Rather than varying 
the /-values, our variable parameters were iV( 12 CO), which scales the entire band, and T ex , 
which determines its shape. Eleven 12 CO bands exhibit similar excitation temperatures with 
averages of T 10 = 6.1 ± 1.0 K and T 2i o = 5.9 ± 0.4 K. Lower T ex were suggested by the fit to 
a' 14, presumably due to effects of higher r. Although /band depends on T ex , the difference in 
/-values between T ex = 4.2 K and T ex = 6.0 K is < 1%, which we ignore here. One additional 
variable fitting parameter was the heliocentric radial velocity. For the main component, the 
average value from nine 12 CO intersystem bands is 14.7 ± 1.0 km s" 1 , where the 1 a error 
is dominated by STIS wavelength calibration uncertainties. This radial velocity is in very 
good agreement with the K I fits mentioned in § 2. 



4.1. Bands Previously Seen Toward ( Oph 

Table 1 lists results of synthetic fits for 12 intersystem bands detected in our spectrum 
of X Per. For the reference /-values we used the MN /j/j/ /-values of the R(0) lines, which 
agree to ± 1% of / ban d-values calculated for T cx = 4.2 K, except for a'14 and a'17, which 
have low- J' crossings and for which the R(0) /-values are off by < 2%. Each ratio of a fitted 
/band-value over the reference MN value was computed from the ratio of each band's fitted 
iV( 12 CO) over the actual iV determined in § 3. Column 5 lists the resulting band /-values 
toward X Per. For comparison, in columns 6—8 we list ratios of /band from this paper, from 
F94, and from theoretical predictions in Rostas et al. over the values from MN. The Rostas 
et al. values are based on multi-f a calculations while the MN values refer to single-f a results. 
The theoretical ratios in column 8 lack formal error bars, although deviations due to the 
presence of other triplet levels and to absolute scales should amount to a few percent. The 
last column directly compares /-values from X Per to those from Rostas et al. 

The first seven bands listed in Table 1 were previously detected toward £ Oph by F94. 
As discovered by F94, the measured /-value of d7 toward X Per is five times greater than 
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the value listed by MN in their Tables 5 and 14. However, MN list another /-value for dl 
in their Table 7, based on measurements by Eidelsberg et al. (1992); the latter result is 
more in line with both IS /-values and was confirmed in room temperatures experiments by 
Rostas et al. The more rigorous multi-iu calculation reproduces the astronomical /-values 
much better than a single- va prediction because dl is relatively far from its nearest parent, 
A2. Our new /-value is 6% (0.4 a) lower than the prediction by Rostas et al., while the F94 
result differs by 1% only. Consequently, we do not confirm the recent laboratory result by 
Stark et al. (2002) that fband(dl) should be reduced by (29 ± 15)% relative to Rostas et al. 

Our synthesis of the optically thick a'14 toward X Per returns an /-value 6% smaller 
than that listed by MN, versus the F94 result of —24%. The agreement with the newer 
predictions is better: the F94 /-value is lower by 19%, whereas the X Per result is identical 
to that of Rostas et al. According to them, single-f a and multi-f a /-values are also identical: 
8.3 x 10~ 4 . Thus, ironically, the strong a'14 cannot help to distinguish between the two 
methods because its (single) interaction with A4 is extremely predominant. All other five 
bands, a'll, a'17, dl2, e4, and e5, are also found to be in very good agreement with the 
multi-Vyi predictions of Rostas et al., i.e., to within 20%. Of these, the F94 a' bands already 
showed an excellent agreement with MN (< 6%), and are within 7% of Rostas et al. The last 
three F94 bands showed large differences (down to —49%) from MN predictions. All three 
show better agreement with Rostas et al.: dl2 now differs by —16%, e5 differs by —23%, 
while e4 is off by —37% (3.1 a). The F94 /-value for e4 remains their most discrepant result 
(after d7). This weakest band of F94 was affected by an improper continuum placement, 
forcing its exclusion from our derivation of iV( 12 CO) in § 3. Stark et al.'s (2002) result for 
e4, that its /-value is marginally too small by (17 ± 15)% relative to Rostas et al., is very 
similar to the X Per measurement: (—11 ± 15)%. Out of the seven bands under discussion, 
the largest X Per difference relative to Rostas et al.'s predictions belongs to a'17 (+20% = 
+ 1.0 a), a band for which a few weaker J" = 2,3 lines are blended with another intersystem 
band, dl2. 



4.2. Bands Newly Detected Toward X Per 

In our spectrum there are first IS detections of four 12 CO intersystem bands: d8, d9, 
dlO, and el2. Of these, /band was predicted by Rostas et al. for d8, while single-tu /-values 
were listed by MN for d8 and el2. We identified d9 and dlO as intersystem CO features 
thanks to their characteristic shapes. Final confirmation was provided by a comparison with 
laboratory bandhead wavelengths in Table 27 of Tilford & Simmons (1972), or through the 
use of molecular state constants from Tilford & Simmons to compute reasonably accurate 
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line positions. The newly detected bands are shown in Fig. 1. Our wavelength calculations 
with the constants for the d 3 A X state identify the new d—X bands as the F 3 (Q = 1) branch 
of the d 3 Aj state. According to the lists of MN, F 2 (f2 = 2) components are at least ten 
times weaker. 

Whereas d8 was already mentioned by MN, who adopted Eidelsberg et al.'s (1992) room 
temperature /-value, the newer experiments of Rostas et al. led to a refined value. Our 
determination clearly favors the new value of Rostas et al. calculated for an IS temperature, 
since the X Per /-value is —24% (—1.7 a) away from their value, but is —61% (—7.6 a) away 
from the value listed by MN. The strongest of the new bands, d9, has a W\ between those 
of d7 and <il2. Its strong detection contradicts its non-inclusion in the tables of MN, who 
imposed a limiting /-value smaller than inferred here. The third d—X band, <ilO, has the 
smallest secure /-value found in this study, ~6x 1CT 6 . The dlO band also lacks a previously 
published /-value. Toward X Per, it is blended with the red wing of the much wider Si IV 
A1402, which was successfully removed by profile rectification. The new IS detection of el2 
close to A8 confirms to within 7% the /-value given by MN, but it was not included in the 
study of Rostas et al. However, in the last column of Table 1 we compare d9, dlO, and el2 
/-values with unpublished multi-tu calculations (M. Eidelsberg, private communication). 
Two bands agree to ± 12% with these predictions, while the /-value of our faintest band, 
dlO, is —34% (—2.3 a) away from its predicted value. 

The high iV(CO) toward X Per also provides us with a clean detection of a fully resolved 
intersystem band from the 13 CO isotopomer, namely, dl2. This band was suggested to exist 
by MN in the C Oph spectrum of A6 analyzed by Sheffer et al. (1992). For 13 CO, the dl2 
band was predicted to have substantial mixing, borrowing 19% of the A6 /-value according 
to MN, versus the mere 1% borrowed by d!2 from A6 in 12 CO. Since the A8 band of 13 CO 
has r < 1, we could find A( 13 CO) toward X Per directly. Using the four-component model 
from 12 CO, we derived A( 13 CO) = (1.94 ± 0.08) x 10 14 cm" 2 , which yields a 12 CO/ 13 CO 
ratio of 73 ± 12 toward X Per. As can be seen in Table 1, the measured /-value of 13 CO dl2 
is 13% (2 a) larger than the single-f a prediction in MN. Unfortunately, multi-f a calculations 
for 13 CO have yet to be performed. 

5. CONCLUDING REMARKS 

The agreement between the band /-values determined from our measurements of ab- 
sorption along the IS line of sight to X Per and the calculations of Rostas et al. is very 
good. Rostas et al. suggested that multi-i^ computations are superior for intersystem /- 
values whenever the spin-forbidden band is not fully overlapped by the permitted band. We 
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concur, because the 11-band average of column 9 in Table 1 is 0.95 ± 0.16, i.e., the —5% dif- 
ference can be accounted for by known T ex approximations and A— X /-value uncertainties, 
the latter estimated to be 10% by Rostas et al. As a caveat, had we retained the e4 band 
in the calculation of iV( 12 CO) in § 3, the average would have changed to 0.90 ± 0.15 due to 
a systematic adjustment of all /-values by —5%. Individually, nine intersystem bands show 
very good agreement with their predictions: a'll, a'14, d7, and dl2 differ by < 0.4 a (7%), 
while a'17, e4, e5, d9, and el2 agree to within 1 a (20%). Such small differences between 
empirical and predicted /-values may be primarily attributed to assigned observational un- 
certainties (< 17% for the STIS data; 16% for N toward X Per) rather than to the possibility 
that current quantum-mechanical treatments are inadequate. The bands d9, dlO, and el2 
were not included in Rostas et al., but a comparison with unpublished multi-fA calculations 
shows that two /-values agree to within 12%, while dlO, our weakest and deblended band, is 
—34% off. The agreement between empirical results and theoretical predictions assures that 
accurate column densities can be extracted from IS spectra showing very strong absorption 
in the permitted bands of 12 CO. Finally, there is a need for theoretical multi-t^ calculations 
of triplet-singlet mixing coefficients in 13 CO. 

We thank Dr. M. Eidelsberg and Dr. D. E. Welty for providing results before publi- 
cation, and an anonymous referee for helpful comments. The research presented here was 
supported in part by a NASA grant for HST program GO-08622 and NASA grant NAG5- 
11440 to U. Toledo. 
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Fig. 1. — Montage of five newly detected intersystem CO bands toward X Per. STIS data 
are shown as points, while a solid line depicts individual fits for the 4-cloud model of 12 CO. 
Heliocentric radial velocity is shown for the R(0) line of all bands. 
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Table 1. Intersystem Bands of CO Toward X Per 



Isotope 


Band 


W x (mA) 


W* p /W c ° ph 


/(X Per) 




/xp//mn 


/F94//MN 


/roo//mn 


/xp//roo 


12 CO 


a' 11 [A2] 


6.52±0.39 


4.6±0.8 


3.04±0.51(- 


-5) 


1.06±0.18 


1.06±0.18 


1.14 


0.93±0.16 




a' 14 [A4] 


44.67±0.40 


2.6±0.1 


8.31±1.30(- 


-4) 


0.94±0.15 


0.76±0.04 


0.94 


1.00±0.16 




a' 17 [A6] 


6.38±0.29 


6.0±1.2 


3.03±0.49(- 


-5) 


1.19±0.19 


1.01±0.20 


1.00 


1.20±0.20 




dl \A2\ 


4.73±0.21 


5.0±0.7 


1.92±0.31(- 


-5) 


4.78±0.77 


5.00±0.75 


5.07 


0.94±0.15 




dl2 [A6] 


9.62±0.23 


5.3±0.7 


5.26±0.83(- 


-5) 


0.67±0.11 


0.58±0.07 


0.69 


0.97±0.15 




e4 [A2] 


4.33±0.33 


7.2±1.6 


1.74±0.30(- 


-5) 


0.72±0.13 


0.51±0.10 


0.81 


0.89±0.15 




e5 [A3] 


15.99±0.23 


6.8±0.5 


7.70±1.21(- 


-5) 


0.78±0.12 


0.53±0.04 


0.69 


1.13±0.18 




d8 [A3] 


2.09±0.27 




8.63±1.75(- 


-6) 


0.39±0.08 




0.52 


0.76±0.15 




d9 [A4] 


6.89±0.26 




3.42±0.55(- 


-5) 








(0.88±0.14) 




dlO [A5] 


1.30±0.22 




5.65±1.30(- 


-6) 








(0.66±0.15) 




el2 [A8j 


3.09±0.24 




1.52±0.26(- 


-5) 


1.07±0.19 




(0.96) 


(1.12±0.19) 


13 CO 


dl2 [A6] 


4.59±0.19 




1.73±0.11(- 


'3) 


1.13±0.07 









Note. — The nearest "parent" A—X permitted band is given in square brackets. References X Per (or XP), MN, F94, and 
R00 are this paper, Morton & Norcau (1994), Fcdcrman et al. (1994), and Rostas et al. (2000). The /-values are computed for 
Af( 12 CO) = 1.41 X 10 16 cm" 2 , and for Af( 13 CO) = 1.94 X 10 14 cm" 2 . Since R00 did not have /-values for d9, dlO, and el2, the 
comparison is with unpublished results from M. Eidelsberg. 



